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1 Introduction 



Semiclassical quantization can be understood as an interference phenomenon. Wave 
fronts propagate along classical trajectories and build up eigenf unctions if they inter- 
fere constructively. In integrable systems tori form a backbone for all classically allowed 
motions and the conditions for constructive interference yield the well known WKB quan- 
tization rule [I], 0, [J. Furthermore, rather detailed information on wave functions, matrix 
elements and selection rules can be derived||, ||. 

It appears that for chaotic systems periodic orbits play a role similar to that of tori 
in integrable systems 0. Gutzwillers famous stationary phase analysis^] of the trace of 
Green's function provides a link between the quantum spectrum and classical periodic 
orbits. This link has been made explicit for hydrogen in a magnetic field || ||] and 
some molecules [TIJ: a Fourier transform of the spectrum reveals sharp peaks at the 
poriods of classical periodic orbits. However, periodic orbits are much too numerous to 
provide a one-to-one connection between individual paths and quantum eigenvalues]]]]. 
Mathematically, this is reflected in the inherent divergence of the formal Gutzwiller trace 
formula [|TT|, [12[] . Experiments on microwave resonances [|13|] show that such relations hold 
for more general wave phenomena as well. 

Methods to overcome such divergences have been developed in the context of general 
dynamical systems fl4], |l5 , |l6|j , where invariant sets can be characterized by their periodic 
points. At the heart of these developments has been the observation that classical periodic 
orbits are strictly organized, both topologically and metrically, and that this organization 
can be exploited to rewrite ill behaved sums over periodic orbits in a convergent form. 
The final result is simple and computationally efficient, sometimes showing faster than 
exponential convergence [|17 1 . 

At present, for a succesful application of the program, a symbolic organization of the 
dynamics is necessary. Because of this requirement, it has been carried out for a few 
systems only, most notably 2-d billiards formed by three disks or four hyperbolas |T8|, 
|T9| , p0| , the anisotropic Kepler problem[|l^, |21| and collinear Helium|22|]. These studies 
(many relevant contributions are collected in a recent conference proceeding|23|) have 
demonstrated that indeed the trace formula does yield semiclassical approximations to 
the eigenvalues. They have also shown that the convergence behaviour can be improved 
dramatically by imposing a functional equation [|T^, P0| . The existence of this functional 
equation is suggested by properties of an S'-matrix approach to quantization [^4], |25], |26 



and by the existence of a similar relation for the Riemann zeta function^, [28, |29| and for 
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Selberg zeta functions on surfaces of constant negative curvature [j30| , ^TJ, but it has not 
been possible to derive it within the semiclassical approximation (compare the discussion 
in Ref. ||) 

Another point discussed in these notes concerns the extension of Gutzwiller's theory 
to include matrix elements as well[33, In principle, this requires the full Green's 
function, which can be expressed as a sum over recurrent orbits, i.e., orbits returning 
to their initial positions in projections, not necessarily in full phase space[^5|. It turns 
outpS] that in the case of sufficiently smooth operators, one can again arrive at ex- 
pressions involving periodic orbits only. This then establishes a complete link between 
experimentally accessible spectra and classical periodic orbits. 

The key technical step will be to express sums over periodic orbits as products over 
periodic orbits. Such products are termed 'dynamical zeta functions' |37|], in analogy to 
Riemann's zeta function p8, 39], which can be written as an infinite product over prime 
numbers, 

primes 

Dynamical zeta functions look similar, the product extending over contributions from 
primitive periodic orbits (labelled p), 

GUs) = E(l-t p (s)). (2) 

p 

The fact is, though, that the objects one is interested in are the products as given 
above, almost never the inverses thereof. Thus there is an inverse relationship in the 
behaviour of the two types of functions. For instance, in accord with one approach to the 
Riemann hypothesis [|39| 1, one would like to identify the zeros of (r(s) along the critical line 
s — 1/2 + it with eigenvalues of a quantum system. However, dynamical zeta functions 
Cd{s) have poles rather than zeros at the positions of the eigenvalues. 

The theory will be discussed for two degree of freedom systems only, since they are the 
at present most interesting class. In many cases, extensions to more degrees of freedom 
are possible. 

The organization of the paper is as follows. In the next section, I will summarize 
Gutzwiller's theory for the spectrum of eigenenergies and extend it to diagonal matrix 
elements as well. The derivation of the associated zeta function is given (2.2) and the 
identification of suitable scaling variables discussed (2.3). 

In section 3 tools necessary for the organization of chaos will be discussed: symbolic 
dynamics (3.1), the connectivity matrix (3.3), the topological zeta function (3.4) and 
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general transfer matrices and zeta functions (3.5). Although illustrated for the case of 
hard collisions in a billiard, the symbolic dynamics can be extended to 'smooth collisions' 
in smooth potentials (3.2). 

In systems with discrete symmetries, zeta functions factorize into zeta functions on 
invariant subspaces. This symmetry factorization and the associated reduction in sym- 
bolics is discussed in section 4. 

The ideas developed here are illustrated for the example of a free particle reflected 
elastically off three disks in section 5. Methods to find periodic orbits (5.1), the conver- 
gence of the trace formula (5.2), the semiclassical computation of scattering resonances 
(5.3), the convergence of the cycle expansion (5.4) and methods to obtain eigenvalues of 
the bounded billiard (5.5) are discussed. 

The relevant parts of a classical periodic orbit theory are developed in section 6.1, 
including a discussion of escape rates and the Hannay-Ozorio de Almeida sum rule (6.2). 

Finally, the issue of semiclassical matrix elements is taken up again and applications 
to experiments are discussed. 



2 Semiclassical periodic orbit theory 

Following Gutzwiller|| 0, a connection between periodic orbits and quantum properties 
can be derived from a stationary phase evaluation of Feynman's path integral. Usually, 
only the trace of Green's function is evaluated, but as will be shown below, a simple 
extension allows for the computation of matrix elements as well. 



2.1 Expressions for trGA 

Starting point is the quantum mechanical expression for the trace of Green's function 
times some observable, Qa{E) = trGA. Expanded in the (complete) energy eigenbasis 
(states \n) of energy E n ), this expression takes on the form 

/ \ , (n\A\n) , . 

so that 

p A (E) = --lm g A (E) = £>|A|n) 6(E - E n ) . (4) 
Thus pa has poles at the quantum eigenvalues, with residues given by the matrix elements. 
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The way to obtain the semiclassical expression for (Q) is to consider a semiclassical 
approximation to the propagator, to Fourier transform to find Green's function and then 
to take the trace. Technically, since the observable can also depend on momenta, one has 



to use a phase space representation such as Wigner's function [40, 36 



Contributions to the trace come from two sources: from the very short paths, where 



the propagator turns into a delta functional], and from the periodic paths. To evaluate 



the first part, one uses a Taylor series expansion of the trajectory in powers of time t, 
and exploits the smallness of t in evaluating integrals 4~3|. This then gives a smoothly 
varying contribution to 

9A,o = J ^S(E - H(p, q))A(p, q) , (5) 

i.e. the average of the observable over the energy shell. N is the number of degrees of 
freedom. Higher order corrections, similar to the boundary and curvature corrections to 
Weyl's law for billiards, can also be calculated^, 

The second part is obtained by approximating the propagator as a superposition of 
contributions from all paths. After a stationary phase evaluation of the Fourier transform 
one obtains the oscillatory part of Green's function[[|T|], 

- ffi(2 , ffi ;„- 1)/2 £i^i' /2 ^' g,/ °-^ /2 . «» 

where the sum extends over classical paths p connecting qx and q2 at a fixed energy E, 
irrespective of the time it takes; S p is the classical action J p c?q, 

/ d 2 S p d 2 S p N 



Ds = det ft* |g (7) 

^ dq'dE dE 2 > 

is the determinant of second derivatives and the index v' p counts the number of caustics 
on the energy shell. 

The next step is to take the trace of GA, 

g A = J d N qG(q,q-E)A(q), (8) 

where for simplicity an observable depending on positions only has been substituted. 
In spirit with the semiclassical nature of the entire calculation one would also like to 
evaluate this integral in stationary phase. This is possible, if A varies slowly on the scale 



of a wavenumber. The calculation then continues very much as in Gutzwiller's case[0, p6 

The phase is stationary if the final and initial momenta coincide, which is the condition 
that the trajectory be periodic. In the neighbourhood of every closed path a coordinate 
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system with q\ along the path and 92, ■ ■ ■ , 9at perpendicular to it may be introduced. 
Using the factorization of the determinant D$ and the fact that up to second order in 
the deviations from the trajectory the action only depends on the stability matrix of the 
classical path, one finds 

d<h---da N \D* |i/2 e tfj.(q)/R-fcW2 = _± e _ (g\ 

{2mh)W-W J Q2 QN 1 Spl \ gi \ | det(M p - 1)|V2- ^ 

where S p is the action along the periodic orbit, M p is the stability matrix around the 
orbit and the phase shift fi p is the Maslov index of the periodic orbit |45[ . 

Since the stability matrix is independent of the position along the path, there remains 
the integral / c?g 1 A(q)/g 1 , which by dq/q = dt may be written as a time integral over 
one period. Allowing for multiple traversals of a periodic orbit, we finally find for the 
contribution of one periodic orbit to gA,osc, 

^ oo g(iS p /h—i/j, p TT /2)r 

9A, P = yA>E | det(M r _ -Qll/2 » ( 10 ) 

with A p the integral of A along the orbit. 

Combining the smooth part (j|) with the contributions from all periodic paths, one 
finds 

Pa{E) = --Imti g A (E) = p A , {E) + J2pa, p ( E ) 

71 

p 

dqdp 



I^A( q ,p)5(E-H( q ,p)) 



+im 7r^^|det(M;-l)|V2 e W 

where 

A p = I P dtA(ci P (t),p p (t)). (12) 



In the final formula general observables A(p, q) have been admitted. The momentum 
p(t) is then fixed to be the momentum along the path at q(t). The average of the observ- 
able over one period of the classical trajectory (p(£), q(t)) is A p /T p . The key requirement, 
beyond the applicability of a semiclassical approximation, is that the observable be suf- 
ficiently smooth. Otherwise a steepest descent approximation to all integrals has to be 
used. ^From the above discussion (eq (|||)) one expects this expression to show poles at 
the (semiclassical) eigenvalues, the residues being the matrix elements. 
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2.2 Selberg's and other zeta functions 



For the density of states the operator A = 1 and thus A p = T p . Then the contributions 
from periodic orbits to (O) may be rewritten as the logarithmic derivative of an infinite 



product over periodic orbits |46|], similar to the Selberg zeta function ||47|| in the theory of 
geodesic motion on surfaces of constant negative curvature . 



With A = 1 and A p = T p the period, the contribution from periodic orbits to (|TTj) 
can be written 



T 



Ra ' osc nh ^ § | det(M; - 1)|V3' 



j(iS p /h-i)j,pTr/2)r 



(13) 



The degrees of freedom enter in the size of the linearization perpendicular to the orbit. 
For a two degree of freedom system, M p is a 2 x 2 matrix of determinant one. If the orbit 
is unstable, the eigenvalues are A p and 1/A P . The denominator can then be expanded in 
a geometric series PHI, 



det (MI 



r 



-1/2 



|A p |-/ 2 (l-l/A p )^ = 5:|A p |- r / 2 A7 r 

j=0 



so that 



| oo oo 

71 " p j=0r=l 



(14) 



(15) 



Upon summing on r one finds 

Ra,osc 



1 OO ,(j) 



p j=o 



l-t 



(f) 



(16) 



where & = e lSp ^ n l ^ p7T ^ 2 \A p \ 1 / 2 A j . Using the relation T p = 8S p /dE, one can write the 



(17) 



quotient in ( |i6| ) ClS db logarithmic derivative, 

Ra,os 
so that finally 



1 00 f) 



Ra,osc = --^-logZ(E) 
71 oh 



with the Selberg zeta function K7I, HB 



Z(E) = II IK 1 " e lSp/h - lflp7T/2 \A p \- 1/2 A- j ). 

j=0 P 



(19) 
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Deriving zeta functions for the other traces involving matrix elements requires a little 
trick and will be postponed until later (section 7.1). 

Depending on ones application, it sometimes is convenient to think of (|19D as an 
infinite product of dynamical zeta functions, Z = JJj 1/ Q with 

VCj = Hi 1 - e^-^^IApl^A^') . (20) 



The leading order term l/(o is Gutzwillers original approximation |7], 48], obtained by 



replacing the determinant in the denominator by the dominant eigenvalue. The discussion 
of convergence (section 5.2) will reveal that indeed this first term is the dominant one. 
All other zeta functions converge absolutely and therefore cannot give eigenvalues or 
resonances near the real energy axis [f49l . 



2.3 Scaling variables 



In general, the actions of periodic orbits are complicated functions of energy [5C|. If the 
Hamiltonian describing the system has scaling properties, e.g. if it is a sum of squares of 
the momenta plus a homogeneous potential, V^(Ax) = X K V(x), then by a virial theorem, 
the action scales with energy like 

s(E)= ^ T(Eo) iE.r^\ (21) 



(2 + «) v UJ \E , 

where Eq is some reference energy and T the period. It thus becomes linear in the variable 

/ J? n (2+k)/2k) 

*=hr) • (22) 



\E / 

In the limit of a billard, k — > oo, A; is essentially the usual wavenumber. Because of the 
simple linear scaling of actions with k, one can use a Fourier transform in this variable 
to uncover the periodic orbit structures |5l|, |J. 



Such scaling Hamiltonians are exceptional. However, in the limit of small H one can 
expand the action to first order in energy around a reference energy Eq, viz. S(E) = 
S(Eq) + T(Eq)(E — Eq). In this limit it is possible to identify an energy interval which is 
classically small (the properties of periodic orbits change little) but semiclassically large 
(the interval contains many quantum eigenvalues). Then approximately E itself is a good 
scaling variable. Formally, this is equivalent to consider the eigenvalues as a function of 
h for fixed classical energy, as used e.g. in the derivation of the spectral statistics of 



integrable systems 52 . 
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3 Organizing chaos 



In this section I provide the necessary formal background on symbolic dynamics, transfer 
matrices and cycle expansions. These tools are important in developing the theory of 
zeta functions and their cycle expansions. Eventually, one might hope to overcome these 
limitations. 



3.1 Symbolic dynamics 



The paradigmatic example of randomness is a coin toss | p3[| , which, at least in principle, 
yields as its outcome a string of heads and tails, with no correlations between consecutive 
events. Thus all strings are possible and equally likely. In a chaotic dynamical system, 
one can find similar behaviour |5^, |55| , p^ , [56] : first dynamics is reduced to a discrete 
map using a Poincare surface of section. Then certain regions in this section are assigned 
'heads' and 'tails'. Depending on where a trajectory crosses the surface of section, it 
will map out a string of heads and tails, and different trajectories will map out different 
strings. As discussed in MacKay's lecture, such a construction is generically possible 



in the vicinity of a homoclinic crossing |)7], [58fl . Different from the ideal coin tossing 
experiment, the dynamical coin is loaded: the probability of occurence of a given symbol 
is determined by the dynamics and need not be the same for all symbols. 

For a certain class of systems it seems possible to extend this symbolic dynamics to 
all relevant regions of phase space. Specifically, for three or four disks arranged in a 
plane so that all lines connecting any two disks are possible (and not shaded by a third 
disk), one has a unique labelling of periodic orbits by disk visitation sequences fl59l, |60 



(see Fig. [T]). As the disks are moved closer together to form a bounded system, orbits 
disappear because of shading by one of the disks|H]]. Nevertheless, one still seems to be 
able to label all trapped periodic orbits uniquely by a string of symbols. 

Here we focus on the three disk system, which is somewhat simpler than the four disk 



billiard relevant for hydrogen in a magnetic field|45[]. In both systems, every trajectory 
can be labelled by the disk visitation sequence. The set of labels assigned to the disks is 
called the alphabet (here: {1, 2, 3}), any string formed from them a word (the trajectory 
shown in Fig. [l] could be labelled by the word 1231312). Evidently, a particle cannot 
bounce off the same disk twice, so that repetitions of the same symbol are prohibited. 
This exclusion of • • ■ 11 • • ■, ■ ■ • 22 ■ ■ • and • • ■ 33 • • ■ is a typical example of a grammer rule. 
The way in which infinite sequences specify periodic orbits is reminiscent of the same 
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construction in the horseshoe map |57], [58|] . In a typical scattering experiment, the ingoing 
direction is fixed and the impact parameter varied. Then there will be an entire interval 
of impact parameters containing trajectories with the same collision future for the next 
n collisions. In Fig. § the interval in impact parameter leading to collisions with disk 
1 is indicated. If the next collision is specified as well, a subinterval is selected. With 
increasing number of collisions, these intervals shrink to a point: thus, there will be 
exactly one impact parameter with the prescribed collision sequence. The past of the 
trajectory depends on the ingoing direction. Repeating then the same procedure for the 
angle rather than the impact parameter a unique value of both impact parameter and 
angle will be identified. This construction is very similar to the way in which strings and 
orbits are associated in the Baker's mapp^ . 



3.2 Smooth collisions 

The previous discussion might seem confined to billiard models. However, it should be 
clear that smoothing the discontinuity at the boundaries of the disks a little bit will 
not change the topology of short orbits. For instance, the potential V(x, y) = (xy) 2 ^ d is 
equivalent to a billiard bounded by the hyperbola xy = 1 for d = and changes to the 
quartic oscillator x 2 y 2 for d = 1. Dahlqvist and Russberg|63|, |64| have followed periodic 



orbits from d = (where a code is known) to d — 1 to establish a symbolic coding for the 
above quartic oscillator. Since the potential becomes more repulsive as d approaches one, 
it is difficult to imagine that new orbits are born rather than existing ones destroyed. 
However, this has to be checked case by case and is not always obvious. 

An alternative way has been proposed in Ref. ||5|, |S5|] . A characteristic feature of 
collisions is that a change in orientation in a local coordinate system takes place. One can 
think of defining the coordinate system using two neighbouring trajectories with parallel 
velocities. During the collision, they will cross in position space, causing a change in 
orientation of the local coordinate system. 

Related to this change in orientation are self conjugate points, where neighbouring 
trajectories started with momentum slightly different from the reference trajectory return 
to it. Then an off diagonal matrix element of the monodromy matrix after a full period 
T vanishes, 

(S*,\ /™_ °V^W (23) 

\bp±J \m xp m pp J \dp ± J 

Regions where one would identify a bounce are bounded by two such conjugate points. As 
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demonstrated in Fig. ^ this method also works for orbits which are very close in position 
space and where it is not immediately obvious whether they undergo a collision when 
approaching the equienergy contour or not. 

For the computation of self conjugate points, one can use the linearized equations of 
motion and the monodromy matrix M. For a two degree of freedom system, this is a 
4x4 matrix, which can be obtained by integrating 16 first order differential equations. 
Two eigenvalues of M after a full period are equal to one, due to the fact that both a shift 
along the orbit and a shift out of the energy shell will be preserved. The interesting part 
of M is the 2x2 matrix m describing neighbouring trajectories in a plane perpendicular 
to the orbit but on the energy shell. As worked out in Ref . J65| , it is possible to introduce 
a coordinate system in which the trivial directions are eliminated and in which closed 
equations for the 2x2 submatrix m can be found. 

An additional advantage of a definition of a symbolic code in terms of self conjugate 
points is the close connection to semiclassics. Since the propagator has an amplitude 



proportional to Jdq±(T)/dp±(0), vanishing of the off diagonal element also signals a 
break down of the semicassical approximation, the presence of a caustic and a change of 
the Maslov phase. This close connection between Maslov indices and the symbolic code 
is important for the cancellation of terms in the cycle expansion (section 5.4). 



3.3 Connectivity matrix 

Given the division of phase space into cells labelled by some alphabet the dynamics 
enters in form of transitions between different cells. The connectivity matrix encodes the 
information whether it is possible to go from one cell to another or not. In its simplest 
form, it is defined by 

{1 if transition from j to i is possible 
if it is not possible 

where the indices are letters from the alphabet. 

For many applications, especially transfer matrices, a generalization based on refine- 
ments of the cells is required. Cells are subdivided and labelled according to the past 
of trajectories in the subcells. All trajectories which cross the surface of section in a 
region inside cell ii, sharing a common past of crossings at i 2 , i 3 , - ■ ■ ,i N , define a unique 
smaller cell which will be labelled ■•%]<!■ When iterated once, they will cross the 
surface in any one of the regions with the label ioV2 ; • • 'i>N-i, where i$ can (in princi- 
ple) be any symbol from the alphabet: the very last symbol is dropped and a new one 
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(24) 



added. The connectivity matrix generalizes to one indexed by words I = i\%2 ■ • -ijy and 
J — 3x32 ■ • • 3n- Since in one iteration only one symbol is dropped and replaced by a new 
one, the connectivity matrix can have entries at positions with coinciding intermediates 
only, i.e. T ItJ ^ for I = ai\ ■ ■ -in-i and J = %\ • • -iN-ib only. 

The number N n of allowed strings of length n is given by the trace of the n-th power 
of the connectivity matrix, 



N n = tvT 



(25) 



For example, for the case of a complete binary code (symbols and 1, no grammer rule), 
the transfer matrix is 



T, 



1 

/l 1\ , (26) 

where the numbers outside the matrix indicate the symbols. One can also define the 
larger connectivity matrix, acting on pairs of symbols, 



bin 



T, 



bin 





00 


01 


10 


11 


00 


f 1 


1 





o \ 


01 








1 


1 


10 


1 


1 








11 


V o 





1 


1 / 



(27) 



It is easily checked that both matrices yield the same strings. The number of strings of 
length n that can be formed is N n = 2™. For example, at length 2, the four strings 00, 
01, 10 and 11 are possible. Evidently, when periodically continued, 00 and 11 correspond 
to the fixed points and 1 of length 1 and 01 and 10 describe the same periodic string. 
The primitive period n p of a periodic string is the length of the shortest block from which 
it can be obtained. By cyclic permutation, there are n p such blocks. 

Thus the total number of strings of length n can be decomposed into the number Md 
of primitive strings of length d dividing n, 



N n = Y,M d . 

d\n 



(2* 



By Mobius inversion ||66|| , one finds 

d\n 



(29) 
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where the Mobius function is defined by = 1, fi(n) = if n contains the square of 
a prime and fi(n) = (— l) k if n contains k prime factors. Some examples are given in 
table 1. 

In case of the three disks with their exclusion rule, the associated connectivity matrix 
reads 



1 

T 3 - d = 2 
3 



1 


2 3 


(° 


1 1\ 


1 


1 




i o) 


gs 


is N n = 



(30) 



with the one for the binary code, except for n = 1 and n = 2 (see table 1). 
3.4 Topological zeta function 

The number of strings of length n is given by trT n and thus dominated by the largest 
eigenvalue of T. The inverse of the largest eigenvalue is a zero of det(l — zT). Using the 
identity det A = exp tr In A and expanding the logarithm, one finds 

/ oo n \ 

det(l - zT) = exp ( - £j — trT n J . (31) 

If T is the connectivity matrix, then trT" = N n . Using the decomposition one can 
then replace the sum on n by one on all primitive periodic orbits p of symbol length n p 
and their repetitions r, 

oo oo z n p r 

= +^ln(l-^), (32) 
p 

so that 

det(l-^T) = H(l-z n v). (33) 

p 

Such products over periodic orbits, formed in analogy to the Riemann zeta function, are 
called dynamical zeta function ||37||, or, if they are derived from the connectivity matrix, 



topological zeta functions. As explained in the introduction they are denoted by 
although it is exactly the product ([33|) and not its inverse which is studied. 



13 



For a complete code on m symbols, the left hand side is easily evaluated to be 1 — mz. 
Note the tremendous cancellations this must imply for the infinite product on the right 
hand side when expanded as a power series in z\ For example, for a binary code, one has 



1/Cwn = 1 - 2z = (1 - zY(l - z 2 )(l - z'Yil - Z y ■ ■ ■ . (34) 

This calculation can be used for the three disks as well. The number of periodic orbits 
agrees with that for a complete binary coding except for n — 1, where there is no periodic 
orbit and n = 2 where there are three rather than just one. Therefore, the topological 
zeta function for three disks can be obtained from that for the binary case, 

(1-z 2 ) 2 



l/Cs-d = 1/0 



bin' 



l-zf 

\2 1 o„2 



= (l-2z)(l + ,z) 2 = l-3z 2 -2z d (35) 

which is still a finite polynominal. If however just one orbit is missing, say one of the 
fixed points, then the topological zeta function is no longer polynominal, 

VCpr = VC&mT^- = 1-Z-Z 2 -Z 3 ---- (36) 

Obviously, the leading zero and thus the topological entropy is still two. 



3.5 Transfer matrices and cycle expansion 

Transfer matrices have the same structure as connectivity matrices and the same van- 
ishing elements, but the l's are replaced by quantities multiplicative along trajectories. 
They provide the connection to classical statistical mechanics |>5], |55], [J7|, |56| and have 
been used in semiclassical mechanics first by Gutzwiller|| |67], Q in his analysis of the 
anisotropic Kepler problem. They also figure prominently in Bogomolny'sf^, [25[] theory 
of semiclassical quantization. Here, the precise form of oj^diagonal matrix elements is 
not so important, since all relevant quantities will be expressed in terms of traces of pow- 
ers of T, which involve periodic orbits only. It should be noted that while there is some 
ambiguity in assigning matrix elements of T, there is none for periodic orbits and thus 
traces of T n : actions, periods and stability exponents are representation independent. 

Entries of the transfer matrix or powers thereof are labelled by the code of the initial 
and final cells. Diagonal elements are thus associated with contributions from trajectories 
that start in one cell and return. An increased resolution with its longer code for the cells 
means that initial conditions for trajectories returning to that cell have to be specified 
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more precisely, collapsing to a point in the limit of infinite resolution. A primitive cycle 
of length n will contribute n times. If p G (d) denotes the different primitive cycles of 
length d and t p the contribution from cycle p, then 

tiT n = Y,J2 dt p- ( 37 ) 

d\n pS(rf) 

Upon substitution into (^Tj) and manipulations similar to the ones that lead to ( PB| ) one 
finds 

det(l-*T) = ~[[(l-z n n p ). (38) 
p 

In general, the transfer matrix will depend on variables (such as the wavenumber in case 
of the Gutzwiller trace formula) and one is interested in the zeros of ( j38|) as a function of 
this variable. Therefore, z — 1. However, z is extremely valuable as an auxiliary variable 
when organizing the product. It is only for the final calculation that one puts z — 1. 

The periodic orbits may conveniently be labelled by their symbolic codes. For the 
case of a complete binary code one thus finds 

det(l-zT) = (1 - zt )(l - zh)(l - z 2 t Q x)(l - z 3 t 001 )(l - Am) 

(1 — z 4 £oooi)(l — z 4 ^ooii)(l — ^ 4 ^om)(l — ^ooooi)(l — z ^oooii) 
(1 - A)oioi)(l - ^ 5 toom)(l - zH 0WU )(l - zHoun) .... (39) 

The cycle expansion is now obtained by factoring out the products and arranging terms 
in a power series in z, just as in case of the topological zeta function, 

1/C = 1-2*0-2*1 -* 2 [(*01 -tlto)] 

— z3 [(tooi — *oi*o) — (*on — *oA)] 
— ^ 4 [(*oooi — *0*00l) + (*om — *011*l) 

+ (*oon — *ooi*i — *o*on + *o*oi*i)] — • • • • (40) 

= J2c n Z n (41) 
n 

The important feature to note is that the contributions to and t\ from the two fixed 
points stand isolated but that all others come in groups. In the limit of t p — > 1 the 
connectivity matrix is recovered and as eq ( PS| ) shows, the cancellations among all higher 
order terms are complete. This is the main use of the topological zeta function in this 
context: it provides a back bone of possible contributions to periodic orbit expressions 
and shows how they are organized. What remains to be checked is that the coefficients c n 
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containing long periodic orbits can be grouped so that cancellations similar to the ones 
for the topological polynominal actually take place. 

Alternatively, one can start from (|31|) and expand directly in a power series in z, 
obtaining expressions reminiscent of a cumulant expansion, 

det(l - zT) = 1 - ztrT - y(trT 2 - (trT) 2 ) - y(trT 3 ) • • • (42) 

For instance, the contributions to trT 2 are (T 2 ) o,oo = ^o, C? 12 ) 11,11 = t\, (T 2 )i ,io = tio 
and (T 2 ) 010 i = t w . Thus, the second term reduces to 

-(tr T 2 - (tr T f) = i(t 2 + t\ + 2t 2 01 - (t + U) 2 ) = t 01 - t h (43) 

in agreement with the cycle expansion (flU). 

For the case of three disks with the ternary alphabet with exclusion rules, the zeta 
function is given by 

1/C = (1-Z 2 t 12 )(l- A 13 )(l-Z 2 t 23 )(l-^123)(1-^132) 
(1 ~ ^ 4 ^1213)(1 — ^ 4 ^1232)(1 — ^ 4 ^1323)(1 — ^ 5 ^12123) ' ' ' 
= 1 — Z 2 ti2 — Z 2 t23 — Z 2 tzi — Z 3 tl23 _ ^ 3 ^132 

— ^ 4 [(^1213 — 1 12^13) + (^1232 — ^12^23) + (^1323 — 1 13^23)] 

-^ 5 [(tl2123-tl2tl23) + '--] • (44) 

Again the terms that stay isolated are exactly the ones indicated by the topological zeta 
function. If the disks are place in a symmetric arrangement then there are relations 
between the orbits and the zeta functions factorize and simplify in a beautiful manner to 
be explained in the next section. 



4 Symmetries 

Many dynamical systems of interest come equipped with symmetries. Continuous symme- 
tries usually give rise to conserved quantities by Noether's theorem. Discrete symmetries 
provide relations between trajectories and can be used to decompose phase space and 
dynamics into irreducible subspaces, just as in the familiar case of a quantum system 
with symmetry where eigenvalues and eigenvectors can be determined for the invariant 
subspaces separately || ^9], [7(], |TJ, [Tlfl . 

Basic to this is the observation that a discrete symmetry can act on an orbit in two 
ways: it can map the set of points making up the orbit into itself or it can map it into 
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a different set which then again is an orbit. In the latter case, the properties (actions, 
periods, stabilities) of the orbit are unchanged, so that some factors in (53) coincide. 

The symmetry group of three circular disks arranged on the vertices of an equilateral 
triangle is C^ v , consisting of the identity e, two rotations C3 and Cf by 27r/3 and 47r/3 
around the center and three reflections C12, C13 and C23 on symmetry lines (see Fig. ^). 
For instance, the rotations map the orbit 12 into 13 and 23. Any one reflection maps 123 



into 132 and so forth. Taking just one representative of every degenerate class of orbits, 



the zeta function (Pj) becomes 

l/C = {i-zH 12 f{i-zH l2Z ) 2 {i-zH l2Vi f 

(1 — £ 5 £l2123) 6 (l — Z 6 tl21213) 6 (l — £ 6 tl21323) 3 • • • 
= 1 — 3^ 2 t\2 — 2z 3 ti23 — 3z 4 (tl213 — ^12) — (^12123 — ^12^123) 

^123) 



-Z 6 (6 tl21213 + 3 ^121323 + ^12 — 9 ^12^1213 ~~ t\ 



— Qz (^1212123 + ^1212313 + ^1213123 + ^12^123 — 3 ^12^12123 — ^123^1213) 

— 3Z 8 (2 ^12121213 + £l2121313 + 2 ^12121323 + 2 ^12123123 

+ 2 ^12123213 + ^12132123 + 3 ti 2 ^1213 + ^12^123 

,2 



— 6ti2tl21213 — 3^12^121323 — 4 tl23^12123 — ^1213) — ' " ' (45) 

A further reduction may be achieved by considering orbits whose trajectories are 
mapped into themselves under a symmetry operation. Then the orbit can actually be 
subdivided into irreducible segments, the full orbit being a combination of several seg- 
ments. Similarly, the plane may be divided into a fundamental domain(Fig. |j) and its 
images under the symmetry operations. One can then define a new code based on the 
group elements needed to map a trajectory back onto the fundamental domain. The 
fundamental operation is a reflection every time the particle hits the boundary of the 
fundamental domain. If this is the only reflection between two collisions with the disk, 
a symbol is assigned, but if two are needed (corresponding to a rotation), then the 
symbol is 1. This new code turns out to be binary without any restrictions. In Table 2 
some orbits and their binary and ternary codes are listed. 

For instance, the orbit 123 is invariant under the rotations C3 and Cf. It can be 
pieced together from three identical segments 12, 23 and 31, mapped into each other by a 



rotation. Under a reflection, this orbit goes over into 321, which is just the time reversed 
orbit and therefore has the same symmetry. Its contribution to the zeta function can 
thus be written 



l-ti 23 ) 2 = (l-t?) 2 , (46) 
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where the new label 1 comes from the fact that the orbit can be mapped back into the 
fundamental domain by a rotation. Similarly, the orbit 12 is invariant under the reflection 
0"i2, so it can be pieced together from two segments 12 and 21. Application of the rotation 
produces two more orbits 23 and 31. Its contribution to the zeta function can thus be 
written 

(l-*i 2 ) 2 = (l-*o) 3 - (47) 

Then there are orbits without any symmetry relations, which have multiplicity 6. Finally, 
there are orbits related by time reversal symmetry but no other geometrical symmetry- 
only one member needs to be computed, which then enters with multiplicity 12. 

The transfer operator is a linear operator and can therefore be decomposed into a 
direct sum of its irreducible representations, implying a factorization of zeta functions 
into products of zeta functions for the irreducible subspaces. An explicit construction 
of the transfer matrix based on the irreducible segments of an orbit is possible, but not 
necessary. Of interest are determinants, for which there is an expression involving traces 
only, 



det(l + M) = l + trM + ^((trM) 2 -trM 



2 

+- f (tr M f - 3 (tr M) (tr M 2 ) + 2trM 3 
6 v 

+ --- + l((trM) d ----) . (48) 

d is the dimension of the representation. Since M is essentially a matrix representation of 
the group element under which the orbit is invariant, its traces are given by the characters 



Xa(s) = tr£) a (g), listed in standard tables^]. In terms of characters, we then have for 
the 1-dimensional representations 

det(l-D a (g)t) = l-x a (g)t, (49) 
and for the 2-dimensional representations 

det(l - D a (g)t) = 1 - xa(s)t + \ (xa(g) 2 - Xa(g 2 )) t 2 . (50) 

Specifically, for the case of three symmetrically arranged disks and C^ v symmetry, one 
has two one- dimensional irreducible representations, symmetric and antisymmetric under 
reflections, denoted A\ and A2, and two degenerate two-dimensional representations of 
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mixed symmetry, denoted E. The contribution of an orbit with symmetry g to the l/£ 
Euler product (ffl) factorizes according to 

det(l- J D(g)t) = (l- X A 1 (g)t)(l-XA 2 (g)t)(l-X^(g)t + XA 2 (g)t 2 ) 2 • (51) 
Using the character table for the C^ v group, 





A 1 A 2 E 


e 

a v 


1 1 2 
1 1 -1 
1 -1 



one finds the following contributions from cycles: 



C3, c| 



(1 
fl 



v) 

3\2 
p) 

2\3 
p) 



A x A 2 E 

;i-^)(i-*p)(i+ h+tjf 

;i-tp)(i + tp)(i + otp-t|) 2 



(52) 



where p denotes the symmetry reduced binary code for the segments of the orbit. 

The outcome of this exercise is that the factorization within the A\ subspace is given 
by that of the binary zeta function fl4T|), and that the one for the A 2 case is similar, 
except that the contributions from orbits with an odd number of O's change sign. More 
interesting is the zeta function for the E subspace, which contains a different pattern of 
terms, 



1/0 



1 + zt x + zHl){l - z%){! + ztvao + z 2 t( 00 )(l - z 2 i\o) 



.2.2 



2.2 



-2.2 



2.2 



1 + Ztiooi + Z ^iooi)(l + 10000 + Z t 



10000 t 



1 + ZtlQIQl + Z 2 t\ Qlm 



)(1 



^ 2 ^ioon) 2 



\ + zt x + z 2 {t\ - tl) + z 3 (t 001 - W 



t-l tr 



+z q 



+z° 



toon + (^ooi ~ 

^00001 + ^01011 



Mo)^i - *oi 



2^00111 + (^ooii — ^01)^1 + (^1 — ^l)t 



100 



+ 



(53) 



Similar decompositions hold for other symmetry groups WQ, 71 . 
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5 The three disk system 



The formal developments of the previous sections will now be applied to the three disk 
billiard ^T], [JI], |75fl . This system is ideally suited to test the methods since a good 

symbolic dynamics is known and tuning of a parameter allows one to study the transition 
from an open strongly chaotic system to a bounded one. 



5.1 Periodic orbits 

The example we will consider is motion of a point particle in the plane with three circles 
removed. The system is characterized by the ratio d/R of the distance d between the 
centers of the disks and their radius R. If the distance between the circles is larger than 
the radius, d/R > 2, then all of the plane is classically accessible and we have a scattering 
geometry. If on the other hand the three disks touch then they enclose a tipped region 
which we will refer to as the bounded billiard. 

The classical dynamics of this system reduces to an exercise in geometry. Several 
coordinate systems are possible: either position along the circumference of the disks and 
parallel momentum or position and length of the segment between any two collisions |60 



or scattering coordinates, i.e. impact parameter and ingoing angle. More important 



is the choice of a numerically stable routine to find the orbit |74], ffql - Drawing on general 
experience in numerical mathematics a multipoint shooting method suggests itself. All 
intermediate point trajectory are allowed to vary, so that for an orbit of symbol length n 
one has 2n variables. A Newton-Raphson iteration will typically converge very rapidly. 



5.2 Convergence of the trace formula 

The exponent in the semiclassical expression for the density of states (|Tl|) becomes for 
billiards S(E)/h = Lk with k = \f2mE/U the wavenumber and L the geometrical length 
of the paths. All lengths can be taken relative to the radius of the disks. The semiclassical 
limit h — > corresponds to k — > oo. 

In many systems, the number of periodic orbits in ([□]) increases and their weight 
decreases exponentially with period. The balance between these effects determines con- 
vergence ]TT], For billiards, the period of an orbit is proportional to its length, so that 
one has equivalent statements for the proliferation of orbits with increasing geometrical 
length. The absolute convergence of the series (|TI|) is determined by the sum over absolute 
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-l/2 e -sL p 


l-\A p 


— 1 /2 g — s Lp 



values. Allowing for complex wavenumbers k = k r + is, this becomes 

OO T 

V V — e~ rLpS (54) 

^^|det(l-M;)|V2 e • ^ 

This expression will converge for sufficiently large s, but will diverge for small s. As we are 
interested in the asymptotic behaviour for long orbits, we can replace det(l— MT) ~ |A p | r , 
where |A p | is the expanding eigenvalue of M p . Then 

oo 

T,Y, L p\ A p\~ r/2e ~ rsLp = E 

p r=l p 

= E^Jo8(l-IM- 1/2 ^) 

= ^logr'oo (55) 

with the zeta function 

VC = II (l - |A p |- 1/2 e-^) . (56) 
p 

This is yet another example of a dynamical zeta function, where the weights assigned to 
an orbit are t p = |A p | -1 / 2 e -sLp . The abscissa of convergence then emerges as a zero of 
this zeta function. 

Fig. |5| shows —So, the negative of the abscissa of absolute convergence. This repre- 
sentation was chosen because of the similarity between — s and the classical escape rate 
7, to be discussed in section 6.2. One curve was calculated using only the contributions 
from the fixed points in a cycle expanded form of (|56), i.e. it shows the zero of 

1 - |A o r 1/2 e- sLo - \A 1 \-y 2 e-' Ll . (57) 

For the second curve the contributions 

- |A 01 |- 1 / 2 e- sZ ' 01 + lAoAxl-^V^"^ (58) 

from the period two orbit have been included. Both curves agree well for d/R > 5. For 
sufficiently separated disks the abscissa of absolute convergence is negative, indicating 
convergence of the trace formula along the real axis. Near d/R ~ 2.8 the curve crosses 
the real axis and approaches ~ 1.5 for the closed billiard d/R — > 2. 

The higher order zeta functions 1 / Q in the Selberg Zeta function contain more powers 
of A J , so that their abscissa of convergence satisfy Sq < Sq ^ < • • • < s . Since none if 
these functions can have zeros in the half plane Imk > Sq , the resonances closest to the 
real axis in the interval < Imk < come from 1/Co alone. If < 0, that is for 
d/R> 2.8 then all resonances are semiclassically bounded away from the real axis. This 
gap has previously been identified by Gasparcl[[73|j. 
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5.3 Scattering resonances 

Scattering processes can be described using the ^-matrix. For the three disk system, 
the explicit expressions of the outgoing waves in terms of the ingoing waves using the 
scattering matrix have been given by Gaspard|73[. Resonances are related to complex 



zeros of tr 5^^§, which is the extension of the density of states to scattering systems [ fTcf] , 
It therefore is given by Gutzwiller's sum over classical periodic orbits, eq ( |TTp with A = 1. 

As discussed in section 4 on symmetries, the Selberg zeta function factorizes into three 
infinite products for the three subspaces of the symmetry group C^ v of the three disks. 
In the symmetry reduced symbolic code, the cycle expansion for the A\ subspace is given 
by the cycle expansion for a complete binary code, eq (f41~p. Using all periodic orbits 
up to symbol length 5, in total 14 orbits, one finds the resonance spectrum (complex 
zeros) shown in Fig. Evidently, the semiclassical calculations agree well with the 
quantum results. Also, allmost all resonances are below the limit [|73| predicted from the 



semiclassical calculations. That the two resonances with smallest real part are above the 
semiclassical limit is one of the deviations to be expected in the deep quantum regime. 

5.4 Convergence of the cycle expansion 

In recent years beautiful arguments for the convergence of the cycle expanded zeta func- 
tion in the realm of semiclassics have been developed. Most important are attempts to 
obtain quantization conditions from a scattering matrix [p4], p5| , 2B|. The basic observa- 



tion is that for finite Planck's constant only a few scattering channels are open so that 
upon neglect of evanescent waves one deals with a finite scattering matrix. Therefore, 
the expansion of det(l — S) in a series in tr S n as in eq ( |4"8| ) terminates. A different 
argument [27, 2£| speculates about the existence of a Riemann Siegel relation and a boot- 



strapping of the contributions from longer orbits. A Riemann Siegel relation also holds 
within the scattering matrix formulation. 

However, from a general point of view, the cycle expanded product should con- 
verge because of the compensations induced by approximations of long orbits from short 
ones[|n], 13, [LlJ. In its most primitive form, the compensation argument applies to terms 
in eq ([|l]) of the form t a k b —t a t a k-i b , involving a long orbit a k b and two approximands a and 



a b. This orbit is a periodic approximation to an orbit homoclinic to a. Substituting 
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the form of t a k b from eq (]20|) for j = 0, one finds 



t a k b -t a t a k-i b = t a k b 1 - e i(S«+S *- lr S aH )/S e -k(Ma+^ fc -i r ^j)/2 



A a A a *:-i 6 


-l/2\ 


A a fcfe 





.(59) 



Because of our definition of the code in terms of selfconjugate points (section 3.2) the 
Maslov indices cancel, /i a + fi a k - 1 b ~ ^a k b — 0- Furthermore, since with increasing k, seg- 
ments of a k b become closer to a, the differences in action and the ratio of the eigenvalues 
converge exponentially with the eigenvalue of the orbit a, 

S a + S a k-i b - S a k b pa const A~ k (60) 
|A a A a fc-i 6 /A afc6 | « exp(-const A- k ) (61) 

Expanding the exponentials one thus finds that this term in the cycle expansion is of the 
order of 

t a H ~ tatak-i b ph const A~ k . (62) 
The number of terms in every order of the cycle expansion is even larger than the number 



of periodic orbits[71|. However, the compensations reduce the size of the contributions 
from the periodic orbits, inducing convergence. In the case of the three disks, compen- 
sations should be best for k = since one then does not have to worry about phases. As 
Fig. ^ shows, compensation is very good for d/R = 6 and d/R = 3. What changes as one 
approaches the closed billiard d/R= 2 is that the convergence is no longer as rapid and 
that pruning and missing orbits cause non monotonic variations (see below). 

The above discussion also shows that with increasing energy more orbits are needed 
to obtain convergence: since the actions increase with energy, so do their differences. 
But if the differences become larger than 0(h), it is no longer permitted to expand the 
exponentials and these terms have to be kept. However, for sufficiently long orbits, the 
compensations will again take place. 



5.5 Eigenvalues for the bounded three disk billiard 

Turning to the bounded billiard, one has to worry about additional problems associ- 
ated with pruning of orbits. As the disks are moved closer, some orbits annihilate and 
dissapear[|6"T|]. For instance the orbits 000001 and 0000011 exist for d/R > 2.016... only. 
Other orbits in the family n l and 0™11 with n > 5 vanish earlier. In addition to these, 
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the orbit = 12 vanishes when the disks touch. These are the only pruning rules for 
orbits of period < 10. 

As the disks move closer together, one would expect the resonances to move up towards 
the real axis, lining up on it for the closed system. The numerical calculations show this 
behaviour only approximately. Including all orbits up to length 10 one finds that the 



resonances lie near the real axis, but not exactly on it^TJ. To remedy this, one can 
speculate about the existence of a functional relation for the Selberg zeta function which 
would put it on the real axis[fM ES . 



Po(^) = -^log(e-^), (63) 



The existence of a functional equation is suggested by the behaviour of the quantum 
mechanical integrated density of states N(E): it is a piecewise constant function that 
jumps by 1 at every eigenvalue. To derive a semiclassical expression for this note that 
the smooth density of states in eq (|TT|) can be written as the logarithmic derivative, 

d_ 
ttOE 
with 

N (E) = J dEp (E) (64) 

the integrated smooth density of states. Together with the Selberg Zeta function (|i~8|), 
( |i9|) one finds 

p(E) = -^\og(e-^Z(E)) (65) 

or for the integrated density of states 

N(E) = -- log (e- i7TNo Z(E)) (66) 

A piecewise constant function can be obtained if the argument of the logarithm is real. 
Then N(E) is constant inbetween zeros and jumps by one, if the argument has a sim- 
ple zero, for then the phase jumps by n. Therefore, semiclassical approximations to 
eigenvalues can be obtained as the zeros of 

D(E) = e~ i7TN °Z(E) . (67) 

This expression is very close to the functional determinant det(E — Ei), (see, e.g., 
Ref. 



For the billiard, the mean integrated density of states is given by the Weyl expansion f42 
containing the area, circumference and curvature of the billiard plus a term from the sin- 
gular tip at the touching point of the disks. The contributions from the periodic orbits to 
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the density of states in the A\ symmetry can be expanded as in eq (f4T|). Fig. |] shows two 
approximations to D(k) vs. k for the three disk billiard, one involving all orbits up to 
length 2 and length 3, respectively. The overall agreement is good and most eigenvalues 
are resolved. Again there are (quantum) deviations for small k. 



6 Classical periodic orbit theory 

When transferred to classical mechanics, the manipulations that lead from the quantum 
propagator to the response function yield an expression for the classical spectrum in 
terms of periodic orbits|78|]. This is useful in interpreting resonances in classical corre- 
lation functions. In addition, these classical expressions contain information on classical 



escape rates in open systems[81] and the Ozorio de Almeida-Hannay sum rule in bounded 



systems [52 



6.1 The classical trace formula and associated zeta function 

Without much difficulty one can carry over Gutzwiller's arguments]^, |7j from quantum 
mechanics to classical mechanics and derive a periodic orbit representation for the clas- 
sical propagator. Starting point is Liouville's equation for phase space densities /, 

f = -iCf, (68) 
with the formal solution 

ft = e- iCt f . (69) 
The equivalent of a Green's function may be defined by Fourier transform in time, 

G d = f dte iut e- iC \ (70) 



and the trace of Green's function then gives the density of states, 

Pd = trG d = J2—^—- (71) 

t 0J — OJa 



The spectrum of an ergodic system |8q, [84]] has one eigenvalue at u = 0, connected with the 
invariant measure, superimposed on a continuum. This continuum will show resonances 
which can be thought of as complex poles uj{. 
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The calculation proceeds from an expression of the classical propagator as an integral 
operator, 

/ t (x) = |dq5(x-q t )/o(q), (72) 

where qt denotes the point q propagated forward in time. Upon taking the trace, the 
delta function singles out orbits that return to their starting point: periodic orbits. Thus 
the delta function only gives a contribution whenever t equals a period T p . Changing to 
a coordinate system along the orbit and perpendicular to it one can do the perpendicular 
integration, resulting in a factor 1/ det(l — M p ), and the integration along the orbit, 
resulting in the period. In view of the delta function the Fourier transform in time is 
easily done, leading to an expression 

oo rp 

^SSsii^^' (73) 

which by manipulations similar to the ones in section 2.2 can be expressed as the loga- 
rithmic derivative of the zeta function 

oo 

Z d {u) = LI IK 1 " \K\- l ^ k e l ^) k+l , (74) 

k=0 P 

where the origin of the power k + 1 may be traced to the fact that in contrast to the 
quantum case one has a full determinant rather than a square root in the denominator. 

6.2 Escape rates and sum rules 

I shall not give a complete discussion of the classical periodic orbit expression and its use 
in determining resonances |79|, |80f| , except for a discussion of the significance of a zero at 
real s. 

Consider first a scattering system such as the arrangement of three non-touching 



disks. One can then ask for the time a particle will spend in the central region [81|. For 
this hyperbolic system the distribution of these trapping times is exponential, and can be 
characterized by a mean trapping time 7 -1 or its inverse, the escape rate 7. One method 
to calculate 7 makes use of trapped periodic orbits. Roughly speaking, if a trajectory 
makes n bounces, then it will be close to some periodic orbit of that period. The orbit is 
exponentially unstable and the probability of staying near it is given by the inverse of the 
expanding eigenvalue. The average over all periodic points, including multiple traversals, 
then is 

00 

EE^IM"^. (75) 
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which is exactly of the form of the semiclassical expression above, with the exception that 
now a full eigenvalue appears rather than its square root. Therefore, by manipuilations 
similar to the ones in eq (|55|), one can write this a the logarithmic derivative of the zeta 
function 

VC = II (l - \^p\' 1 e- T ") • (76) 
p 

Using cycle expansions up to n = 1 and n = 2, respectively, one finds the curves shown in 
Fig. p. The agreement between the two curves is an indication of the rapid convergence 
of the cycle expansion. Also shown are the results of a classical Monte Carlo simulation 
of Gaspard^O . 



There is a direct connection between the escape rate|81j| and the classical sum rule 



proposed by Hannay and Ozorio de Almeida p2|. According to one formulation (see in 



particular section 17.7 of Gutzwiller's book||), the sum 

V Tp ~ i (77) 

T ^ T det(l-M p ) ' 1 > 

so that the set of periodic orbits approximates the phase space measure in the sense 
that every average over phase space can be expressed as an average over periodic orbits. 
Rather than defining a set of orbits by their periods, one can also take all periodic orbits 
of a fixed symbol length. Then the equivalent statement is that 

V ^ « 1 . (78) 

This may be expressed as a logarithmic derivative of a zeta function by multiplication 
with z n and summation over n. That the sum approaches a constant then means that 
the zeta function 

i/c = na - iA P rv*) (79) 

V 

has a zero for z = 1. But this is just the escape rate zeta function eq (|76| ) for s — 0, 
so that a vanishing escape rate indicates that the periodic orbits represent the invariant 
phase space density (in hyperbolic systems). 



7 Matrix elements 

In this final section, I will discuss the derivation of a zeta function for the diagonal matrix 



elements and indicate the extension to off-diagonal matrix elements [SB . 
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7.1 Diagonal matrix elements 

As discussed in section 2.1, the final expression for trGA with some sufficiently smooth 
observable A becomes: 



Pa(E) = J -^-A( q ,p)5(E - H( q ,p)) 



7 00 A 

_|_Ttti W p JiS v (E)/h-iwy/2)r I oq\ 

vr^^|det(M;-l)|V2 e ' ^ 

where 

A p = [ Tp dtA{q p {t),p p {t)). (81) 



The derivation of a zeta function for matrix elements is similar to the calculation in 
section 2.2, except for one step. One first arrives at 

1 00 A fa) 

^^EE/i]). ( 82 ) 

where t^> = e lSp / ft_l ' 1 * ,7r / 2 |A J ,|~ 1 / 2 A~ J '. In the case of the density of states, the observable 
is A = 1 and the prefactor A p = T p = dS/dE so that one can immediately write the 
quotient in (|T6D as a logarithmic derivative and thus arrive at the zeta functions. Here, 
one may proceed similarly after introducing an auxiliary variable 7 and an extended 
t (j) = Since = ^(7 = 0) and Apt® = dt p ^/dy\ 1=0 , one can write 

^.c=^EE^log(l-t«)U (83) 
and thus finally 

i f) 

(84) 



R A , osc =^^Z^E) 



7=0 

with the extended Selberg type product 

00 

Z(y,E) = [] - e -^e^ /?i -^ /2 |A p |- 1/2 A7). (85) 
v 3=0 

For 7 = this dynamical zeta function coincides with the Selberg Zeta function Z S (E) 
for the spectral density p(E), 

Z s (E) = Z( 1 = 0,E). (86) 
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Eq (HH) may again be evaluated using the cycle expansion. Consider the case of a 
hyperbolic system with complete binary symbolic dynamics. Labelling the contributions 
from the periodic orbits by their symbolic code, one has to consider products of the form 



Z( 7 , £) = (!- t )(l " *i)(l " *oi)(l " *ooi)(l - toil) 



(87) 



which expand into 



Z( 7 ,£) = l-to-ti 



(toi — tati) — (tooi — to^oi) — (ton — titoi) ■ • • 



8 



8 




(89) 



The suggestive argument for convergence (section 5.4) is that a long orbit t a b can be 
decomposed into shorter ones (labelled a and b), which shadow it; therefore, the term 
tab — t a tb is small. However, the same can be expected of the observable: if the short 
orbits a and b are close to ab, then A a b should also be close to A a + Ab and the convergence 
properties of the above cycle expansion can be expected to be the same as the ones for 
the density of states. 

7.2 Correlation functions 

A similar result may be obtained for correlation functions. All one has to note is that 
in the semiclassical limit the time evolution of the operator B is given by the classical 
observable B at the time evolved positions. Therefore, one obtains correlation functions 
if the observable A(p, q) is replaced by B(p(t), q(t))C(p, q), where (p(£),q(t)) denotes 
the point reached by a trajectory starting at time t — from (p, q). The final result is 
exactly as above, except that A v is replaced by 



the correlation function of B and C along the periodic orbit. The first term, the contri- 
bution of paths of zero length, becomes the classical correlation function. 
Quantum mechanically, tr G B t C evaluates to 




L 



dt'B(q(t' + t),p(t' + t))C(q(t'),p(t')) 



(90) 



tr GB t C = Y. 



n 




(91) 
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where 



[71 



B t C\n) = ^2(n\B\m) (m\C\n)e i{En - Em)t/n ; 



(92) 



one therefore expects the semiclassical expression to have poles at the eigenenergies, with 
residues the matrix elements (n\B t C\n). 



7.3 Periodic orbit spectroscopy 

The developments in the previous section pave the way for what one may call 'periodic 
orbit spectroscopy' |85|]. A Fourier transform of spectra in terms of a suitable scaling 
variable very often reveals sharply defined peaks, which according to semiclassics can 
be assigned to periodic orbits. The most prominent example of this is hydrogen in a 
magnetic field, where the comparison between classical and quantum mechanics has been 
pushed very far|| |], |86 . 



The formula for autocorrelation functions should be useful in describing Rydberg- 
Rydberg transitions. Calculations by Zeller[^7]] for hydrogen in a magnetic field have in 
fact revealed periodic orbit structures for transitions between highly excited states. A 
quantitative comparison would be highly desirable. 

As an application outside atomic physics, Wilkinson 0] has suggested to apply the 
above formula to conductance problems, since the conductance is related to the Fourier 



transform of a correlation function by the Kubo Greenwood formula |88|, p9|. This relation 
may be useful for the AC behaviour, but must be expected to fail in the DC limit since 



it does not describe correctly the infinite time limit of the correlation functions [36]. 
Periodic orbits also figure prominently in recent interpretations of the spectra of small 



molecules Their presence becomes plausible via the classical periodic orbit formula. 
It would be interesting to identify semiclassical corrections to this. 



8 Final remarks 

Evidently, periodic orbit expansions provide a powerful tool for all kinds of calculations in 
classical and semiclassical mechanics. Almost all quantities of interest can be expressed 
in terms of periodic orbits or, more precisely, traces of transfer matrices. The cycle 
expansion with its arrangement of periodic orbits into compensating groups provides a 
numerically efficient and convergent method to evaluate periodic orbit expressions. 
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The tests presented here and others have aimed at verifying the semiclassical approx- 
imation level by level. This is mainly of theoretical interest, to strengthen confidence 
and explore the range of validity and reliability of the Gutzwiller trace formula, which 
after all is derived using many approximations. The indications are that it is a bona 
fide semiclassical theory, which shows large deviations from exact quantum calculation 
in the quantum domain of small wavenumbers, but better agreement at higher energies 
in the truly semiclassical regime. There still is the question of whether the present form 
of semiclassics can actually give the positions of resonances accurately to better than a 
mean spacing, considering that the propagator includes terms of order h only, whereas 
the mean density of states is 0(h 2 ). Indications are that the error will not be noticable 
until one gets to very high energies. However, it should be more important as the number 
of degrees of freedom increases. 

On the practical side, this appears to be just a minor limitation. The strength of the 
semiclassical expressions is that they can explain long range correlations in the spectra, 
somewhat independent of the fince structure of the spectrum. The aim will be to calculate 
semiclassically from just a few periodic orbits the spectra of small molecules and to use 
this as a tool in analyzing spectra. The most prominent example of this class is hydrogen 
in a magnetic field, but currently some molecules are under investigation where such 
quantitative predictions may be possible as well. 
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M n {2) 
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M n (4) 
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iV(iV 2 - l)/3 
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8 
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N 2 (N 2 - l)/4 
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(iV 5 - N)/5 
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(iV 6 


- N 3 - N 2 + N)/6 
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670 
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(N 7 - N)/7 
N 4 (N 4 -l)/8 
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312 


2340 
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30 


810 


8160 
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iV 3 (iV 6 - l)/9 


56 


2184 


29120 


10 




— iV 5 — iV 2 + N)/W 


99 


5880 


104754 



Table 1. Number of prime cycles for various alphabets and grammars up to length 10. The 
first column gives the cycle length, the second the formula ( p8| ) for the number of prime cycles 
for complete iV-symbol dynamics, columns three through five give the numbers for iV = 2, 3 
and 4. 
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Table 2. C 3v correspondence between the binary labelled fundamental domain prime cycles 
p and the full 3-disk ternary {1,2,3} cycles p, together with the C 3v transformation that maps 
the end point of the p cycle into the irreducible segment of the p cycle. The degeneracy of p 
cycle is m p = 6np/n p . The shortest pair of the fundamental domain prime cycles related by 
time symmetry are the 6-cycles 001011 and 001101. 
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Figure captions 



1. The three disk billiard. The trajectory can be labelled by its disk visitation sequence, 
viz. 1231312. 

2. Construction of impact parameters for long collision sequences. The dashed lines indicate 
the interval in impact parameter leading to collisions with disk 1. Initial conditions in 
the two shaded intervals lead to collisions with 1 and 2 or 1 and 3, respectively. 

3. Two orbits for hydrogen in a magnetic field in semiparabolic coordinates. The energy 
is +0.2 Ry. The self conjugate points are indicated. The four disk code of the two 
trajectories is 142124 and 14121214, respectively. The collisions indicated in bold face 
are not immediately obvious in configuration space but clearly identified by the presence 
of two additional self conjugate points. 

4. Symmetries of the three disk scattering system. Indicated are the three reflections Ui 
across symmetry lines, the two rotations C3 and Cf around the center by ±27r/3 and the 
fundamental domain (shaded). 

5. The abscissa of absolute convergence for the system of three disks as a function of the 
distance/radius ratio d/R. 

6. Some resonances of the S"-matrix in the A\ subspace for d/R = 6. The crosses are 
the exact quantum results, the open circles the result from a semiclassical calculation 
involving all periodic orbits up to symbol length 5. Notice that all semiclassical resonances 
lie below the abscissa of absolute convergence, \mk c = —0.123..., but that the first two 
quantum resonances lie above. 

7. Convergence of the cycle expanded zeta function for the three disk system. The terms c n 



contain all contributions from orbits and pseudo-orbits of length n (see eq (41)) at k = 0. 



8. The functional determinant ( |67| ) for the closed three disk system including orbits of 
symbol length n = 2 and n = 3. The dotter vertical lines indicate the positions of the 
exact quantum eigenvalues. 

9. Classical escape rate 7 as a function of the distance between the disks. The crosses are 
the results of a Monte Carlo simulation by GaspardpT]. 
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